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In this paper, a Brownian motion of order n (n> 2) is defined by a probabilistic approach 
different from Hochberg’s and Mandelbrot’s. This process is constructed from sums 

of independent R!/"-valued random variables (rv) (where R!” = {zeC; z'eR,}). Many 
properties of the real standard Brownian motion are generalized at order n, but in the case 
n> 2, it is interesting to describe the Brownian motion of order n on the o algebra 

®@ [B(R!")] R+ [where B(R!”) is the ø algebra generated by sets of type A (0,h) ={zeC; 
z"e[0;A”[}, (heR* )]. This o algebra is totally different from @ [B(R)] F+. Thus this study shows 
the fractal nature of the Brownian motion of order n, and given invariance scale (self- 
similarity) properties. Then, a stochastic integral and an Itô-Taylor lemma at order n are given 
to allow the representation of the solution of the heat equation of order n by a probabilistic 
average. All these results can be obtained via nonstandard analysis methods 
(infinitesimal time discretization). Finally, one remarks that this process has a.s (almost 
surely) continuous sample paths, infinite variance, and independent increments, whereas the 
fractional Brownian motion of Mandelbrot has a.s continuous sample paths, finite 


variance, and interdependent increments. 


I. INTRODUCTION 


All the shown results come from the properties of the 
complex roots of the unity of order m and the set 
Ri/"={zeC;z"eR , }. We want to study the properties of 
the heat equation of order n (1) by a probabilistic way 
generalizing Itd’s calculus, 


Ou 1 ĝu iai 
ðt (x,t) =al Ox" (x,t), (x,t)e +? 


(1) 
u(x,0)=f(x), (xeR); 
with f a bounded and continuous given (RR) function. 
Following Block,' this problem is essentially linked 
up to the following ODE of order n—1 by using the 
Fourier transform on (1). Then, the fundamental solu- 
tion of (1) verifies (2): 


n—1 


at (2) +(n—D-z-w(z) =0, (2) 


where w(0),w‘ P? (0),... w"? (0) are given real con- 
stants. 

With the help of the notion of generalized pictures 
developed by Block,! one remarks that (1) is invariant on 
R!”, That is right, the operators d/dt and (1/n!)-d"/dx" 
accept, respectively, exp(At) and exp[(!A)'“"-a,-x], 
(o,=exp(2imrk/n); k=0,....2—1) (AEC) as a system of 
characteristic vectors. Then, to follow the basic construc- 
tion of the Brownian motion, it seems necessary to study 
sums of independent R!/”-valued rv with the same law in 
each branch of Ri, the result is C valued, and to have 
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finite characteristics on this limit of rv sums, one will 
often examine the rv sums on the tensorial product o 
algebra constructed from the o algebra B(R!”) gener- 
ated by sets of type: A(0,h) ={zeC;z"e[0,h"[}, (heR* ). 
Then, it reveals the fractal nature in the complex plane of 
the sample paths of this Brownian motion of order n 
(n>2). Furthermore, one obtains a weak central limit 
theorem of order n like Hochberg’s,? but in a true prob- 
abilistic context on the complex plane. Of course, the law 
of the rv sums is not infinitely divisible in the Levy-Doob 
sense.’ Noting down X the obtained limit, one has the 
“strange” properties: 


(WAeC), 
E{exp(A:.%)}=exp(A"/n!) 
and 
E{ | F|? }=+ 0. (3) 


This limit is linked up (via a definite isometry), to the 
cone 


+o 
l (R4)=)(@k)keN; aR}; » a< +0}. 


From that, with a linear cut approximation, one deduces 
the construction of a C-valued process X/,), with station- 
ary, independent, self-similar increments and a.s 
continuous sample paths definite from the o algebra 
a [B(R!/”)] F+, which is an infinite product of copies of 
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B(R'/") (in the compactified Alexandrov’s sense). Xin] 
verifies, as well, the following “strange” properties: 


(Y O<s<t< +o), 
(WAEC), EL |X in (t) — Xin (s) |P}= + 0 


and 


A” 
Efexp(a- [Xn] (t) = TOE a-o]. 
(4) 


Finally, X,,; has several properties generalizing these 
known at order 2 (Brownian scaling, isotropy on the 
branches of R”, no differentiability of the sample con- 
tinuous paths, inversion at order n, Levy’s continuity 
Modulus of paths with bounded nth-order variations, 
etc.). From that, one defines stochastic integrals of the 
following type: 


t 
f F(X in (S)] (aX in (8) <k<n), (5) 


for feL,,(dP x dt) or f having an analytic prolongation to 
C, with “good” properties at infinity. Here (5) is defined 
from a prolongation of an isometry on the cone /,(R,) 
via the strictly definite positive functional E'/"{(-)"}. 
Then, via the Taylor formula developed at order n, one 
obtains an It¢-Taylor lemma at order n for a certain class 
of analytic functions on C, 


a 
È GLO Kim OY dX m]! 


di AX in (9) 


d 
exp( dxo g) -ra| -(F% OD 


(130) (6) 


(it is a random drift in the operational sense). 

It is then easy to deal with the solution of the heat 
equation of order n and his representation via the notion 
of heat or Hermite polynomials of order ,°° and general 
problems (Feynman-Kac problems at order n)”® studied 
directly in Ref. 9 with the help of Trotter’s formula. Fi- 
nally, a study of basic measures implicitly used in this 
paper to define the Brownian motion of order n is given, 
and one links together the Hochberg’s “probability”? and 
that described in this paper. This relation is obtained 
from an extension of the Fourier transform, because the 
support of X7,) is of fractal nature, whereas the support of 
Hochberg’s is real. Many properties of X;,) are totally 
different from those of the Lévy-Mandelbrot process.!° 
One can equally remark that all these results can be 
showed via nonstandard analysis methods following 


Keisler''! with his nonstandard Itô’s calculus or Berger 
and Sloan’? with their nonstandard Hochberg’s calculus. 


Il. N-GAUSSIAN RANDOM VARIABLES 


Let (0,7,P) denote a reference probability space and 
n a natural number (n > 2; the cases n= 1,2 are known). 
Noting down @,=exp(2imk/n) for k=O,....n—1; then 
one has the following results. 

Definition 2.1: One calls Rademacher rv of order n, 
all discrete rv €j} with a uniform distribution in N= {a,; 
k=0,...,.n— 1}, i.e., 


P{E[,)=@x}=1/n, for k=0,1,....2—1. (7) 


Lemma 2.2: Denote by 5,9 the Kronecker’s symbol. 
Then one has 


(Y peN), (Y q=0,....2—-1), E{eP*9=8,0, (8) 


(V AcC), Efexp(A-e;,))}=coshy,)(A) 


+ 00 A"K 


= A akl 


(9) 


where cosh;,) is the hyperbolic cosine function of order 
13-15 
n. 
Proof: It results directly from the properties of the 
complex roots of the unity of order n, because of the 
following calculus: 


n—1 
Bet =n hh E opt) 


k=0 


n—i 
=l. k = 
= ( 2, ole) = 20 


(that is a geometric series), therefore the result. 
Notation 2.3: One notes that 


+o 
Ey,(2) = » Tisha (0) (10) 
[with T(x)=sf%e'2-!-dt (x>0)], the Mittag- 
Leffler function of order n.'&'8 
Properties 2.4: From Refs. 16, 17, and 19, one knows 
that £,,,(—x) (x >0) is a function totally monotonic. 
Definition 2.5: Note down X,, the real positive rv 
defined by his Laplace transform, i.e., 


(Y A>O), Efexp(—A-X,)}=£,,,(—A/K,,), 


(K,=(n!)"”"), (11) 
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and €j) a Rademacher rv of order n independent of X,; 
then, one calls standard n-Gaussian rv, all R'/"-valued rv 
noting that G,,; with a law verifying 


Law Law 
Gin) ~ Etn) Xn (one remarks that: |G,,)| ~ Xn). 


Remark 2.6: The density of the | G;,;|’s law has been 
given before by Pollard’? and Zolotarev.'* It has the fol- 


lowing expression: 
te sin(ka/n) k 
5 eT (145 | 


k=1 


n 
Pa(x) =z Er 


[Kea (13) 
Lemma 2.7: One has the following relation: 


law 
lGin | ~ [A(Z] 7" (14) 


(where Z, is a 1/n-stable real positive rv normalized by 
his Laplace transform), i.e., 


(Y p20), E{exp(—p-Z,)}=exp(—p'””). (15) 


Proof: See Refs. 20 and 18 for the proof (the duality 
relation between Z, and |Gj,)| ). 
Remark 2.8: For example, at order 4, one gets 


Law 
[Gia] ~ (3)7'*: VIGI 1G", (16) 


where G and G’ are two independent real standard 
Gaussian rv; this remark is used implicitly in Ref. 9 for 
the treatment of the Feynman—Kac formula at order 4. 

Lemma 2.9: Here p,, (see Ref. 13) has, equally, the 
following expression: 


n? +0 L1 
Prl x) =n f r exp| =f" 


sin dt. 


2r _ [2m 
—K, cos < )ær K,, sin( 2) 
n n 


(17) 


Proof: This expression is a Hardy function,'® and one 
can see that there exists an inversion relation (called a 
duality relation by Zolotarev'*), between the 1/n-stable 
rv and |G,,,|. To obtain the result, it is sufficient to write 


r(145) = i ented 
n o 


+o . 
=n Í e`" -u*t”-!.du, 
0 


and then to permute the symbols f and = without any 
problems. 

Remark 2.10: This formulation is near Hochberg’s 
density,” which is written for n even as 


1 pa Be 
Pas f, exp) —ixé— -i dé 
1 +o —z" 
=|, a(z) cos(xz)dz. (18) 


Lemma 2.11: Let Gi, be a standard n-Gaussian rv; 
then 


(V AcC), Efexp(AG,,))}=exp(A”/nl). (19) 


Let (GO) pene be a family of independent standard 
n-Gaussian rv and (a) xensx a Sequence of the cone 


+o 


la (R4) = | (œk) kensa > 0; >» a< + 00 |. 


Let S,=24_a,G{,); then S, has the following remark- 
able properties: 


La 
S,~ Sp for k=0,...n—1. (20a) 


(V AeC), Efexp(AS,)} 


(2, e) 


An 
=exp| = Esp). (20b) 


=exp 


But, if 


p 
lim È aj=+o0, 
Ppot+o k=l 


then 


lim E{|S | }= +o% (20c) 


p> +00 


(C is too big a support to correctly analyze the S,’s prop- 
erties; the S,’s support is of a fractal nature). 


Law 
Proof: From (12), one has Gin ~ €[,,X,3 by using 2.5 
and (10), one gets 


(V peN), E{XT?}=(np)!/[ pl(n!)?]. 


Then 
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ELG SH Ele a VEN TH 
=6, 0E 1X7) 


(ip)! 
99 p(n!) 


=ô for all q=0,...,n—1. 


From that, one easily gets (19) by writing the analytic 
expression of exp. By construction (20a) is trivial. Equa- 
tion (20b) results from (19), from the independence of 
the Gi) and from the fact that ea: js =0 for q= 0,...572 
—1. For Eq. (20c) one has E{|S,|"}= = 4 1a; there- 
fore the result. 

Lemma 2.12: In the case when n is even, i.e., n=2p, 
(peN*), one links together Pp and pzp by the following 
“weak” identity. 

For all function g in the real Schwartz class (g has 
an analytic extension to C), 


E{P( Gin) }=Elp(G))} 


—A" 
=— =Í BAdex( SP ) a, (21) 


where P(A) =Sp exp(iAx) p(x)dx and Gi is the Hoch- 
berg “rv,” corresponding to Pye: 

Proof: Since E{exp(AG,,))}=exp(A"/n!), one gets di- 
rectly the result by extension of the Parseval identity on 
the branches of Rin (the idea of the proof is the same as 
for Lemma 2.15). One remarks that the o algebra asso- 
ciated with Gj,) is generated by the sets A, h) and that 
Gin, has a probabilistic sense, whereas Gy is defined on 
B(R), and has a signed density. Thus, one has a weak 
identity compensating the difference of the laws by dif- 
ferent supports. 

Lemma 2.13: Let f denote a bounded and continuous 
(R-R) function, allowing an analytic extension to C, f, 
which verifies 


(A A,M>0), |f(z)|<Mexp(A[z|”) 


(n-'4n'"=1), (22) 


i.e., f is of order n’; then 


lim A7! [E{ F(z+zoh""Gin])} FoI Ie 
ar Z+ Gipi- ed 


More generally, for a holomorphic semigroup exp(zA), 
one gets 


{exp(zo4") }[ f(z) ] = (exp(z}"Gin4) J (z)}, (24) 


and for an analytic function @ in C, 


{exp(zop(A))}[ F(z) ] 
+o p” 1/k 
= (exp ( pA (0% 0) off) 4). 


Fte-+z99(0))) (25) 


(with a formal sense), where (Ge ) cken») IS a Sequence 
of emia tapersie standard k-Gaussian rv (in particular, 
G{i) and Gf) is a standard Gaussian rv), and A is the 
infinitesimal generator of exp(zA). 

Proof (Wh,z,Z9€C): 


A- (ELF (z+ 29h" Gin) )}—F(2)) 
=h7! (2HE FC PEY 


hitin 
eam); 


h 2 
+25 zy ELG in LS (2) +29") 
since E{GfR*%=5,o[(np)!/(pl(n!))], for all peN and 
q=0,...,.2—1, one gets 


xn) 


lim h7"[E{ f(z+29h""G1))} F(z) =z ae 


h=0 


For the generalization to a semigroup, the same ar- 
guments are used. Then, one studies the semigroup 
exp{(t/n!)d"/dx"} from the examination of a central 
limit theorem at order n and a stochastic process gener- 
alizing the standard real Brownian motion. 

Lemma 2.14: Let (ay) ene denote a sequence 
of positive real numbers in the cone /,(R,) 
={ (ak) kns Eak < +o}, (GY) rens a family of in- 
dependent standard n-Gaussian rv and S, the sum de- 
fined by S,=2%_ aG; then, w-limp— + S, exists in 
the following sense: lim, _., E{Q(S, )} exis for all 
functions p, verifying the following conditons: gp is in 
L,(R,dx), continuous, with an analytic extension @ to C, 
such that 


+o 
(v O<x< 2 at), 


=1 


P(A)exp{(—iA)"x/n!}eL,(R,dA) (AeR). 


Proof: The set of g functions is denoted by H,.(C), and in 
all the cases one can choose @ of order n'=n/(n—~—1). 
This kind of convergence is an extension of the usual 
weak convergence.”!” For the proof, one comes down to 
the Stone-Weierstrass theorem of approximation of an 
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analytic function on C by a sequence of polynomials, 
(P;) n in the sense of the uniform convergence on a com- 
pact set of C. 

Let (O<q<p), (p,gEéN*), and (aeR*). Let B, 
={zeC;|z|<a}, and let T, be the stopping time defined 
by 7,={infk>0, (keN*): |[S,|>a}; then one gets 
|Spar,| < @ and since lim, ,,. Ta= + œ, it results that 
|E{P(S, )-9(S,)}|<M (because vpeN*, lim,_. +. o SPATa 
= S,), and we “also come down to a “compact problem” 
on Spar, Furthermore, P(A exp —iA)"x/nl} 

eL,(R,dA), and then, noting that - ŞS, =S,AT, (3 Mp 
M >0), such that IELE, )— ~ GS, \}]<M,<M< +o; 
also, one gets 


(Y e>0) (3l>0), 


|E{@S,) —G(S,)}| <2e+ | EIPS) -PS | 


[because of the uniform convergence of P; to @, 
(J++ œ), on the compact set of C, B,]. 
Furthermore, for j= 1,...,.2—1, E{(S, _§ t= 0 and 


p-1i 
| E{(S,—Sq)"}|< 2 ak +0 
=q 
(q>+o) ((a%)enee/,(R, )). 
One works on a compact set (because of the stopping 


time T,); for all peH,(C) (3 M,>9), |E{p(S,) | 
<Mg< +0; then 


d” 
(3 K,,K>0),Sup Pz) | <K,<K < + œ. 
|z| <a dz" 
Hence 
(Y ¢e'>0) (3 géN), 


|E{@(S,) —G(S,)} | <2e+ Ke’, 


by using the Taylor’s development of P; at order n from 
Sy One then deduces that lim,- 4 «o ELSIS, )} exists for 
all a< +œ, and going to the “imit a— -+ œ, one finally 
gets the result, i.e., lim, + E{p(S,)} exists. 

Lemma 2.15: Under the same hypothesis as in 
Lemma 2.14, one gets 


lim E(G(S,)} 


pote 


-5 tome 


us [ lim E{S"}]}dA. 
pr+o 


(26) 


Proof: Let Sam = (1/nz, -= “oR, ) and p,,(x) be the 
density defined in (13). Then” (Y zeC), G(z) 
= (1/21) SpG(A)exp( —iAz)dA, and one deduces that 


E{p(S,)} 


GERI GEL.) 


P 
xp ( » co Pa "py (Xp)dx,°>-dxp 


1 
R!” Rv" 2a 
p 
xexp ~ia( > ain) aa 
k=1 
Pal ¥1)dx1**“Pn(%p)dXp 


By permuting S RY" and fpr, one remarks that 


Son Spiorwt dr 
R n 


K 
= E{exp(— iA, )} 


=exp{(—ia,)"/n!}. 
Then 


AN $ allan 


k=1 


1 
ESD f, Pep 


Therefore the result, going to the limit (p— + œ ), since, 
E{S,}=Ē2k-1@% In fact, one has used a kind of convo- 
lution on the branches of R”, the Parseval identity, and 
the following result on S,: 


(V AeC), (Y peN*), 


wil 2.) 


Remark 2.16: It is easy to improve this result with a 
greater class of functions g, verifying only (i) pọ is a 
continuous function in L,(R,dx); (ii) 
(Y O<x<dfp%a%, @GlA)exp{(—iA)"x/n!}eL,(R,dA); 
and (iii) (3 6&>0), such that allows an analytic ex- 
tension @ to B(8) ={zeC;|z| <5}. These conditions give 
a sense to E{@(S PAT, )} (o <a<ô). 

Lemma 2.17: In fact, we have 


Ef{exp(AS,)}=exp 
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2 ofA (—iA)"E{S;} 
EIPS) =z f Pe Ae an 


=E{Q( EM "{S5}Gin\)} 


1 n—1 
= lim X E{P ESH Gim DI; 
q= + 00 k=0 


(27) 


The notations and conditions are those of Lemmas 2.14 
and 2.15. 

Proof: For the two first identities, the proof is given in 
Lemma 2.15. For the third, (P q)geN is a sequence of poly- 
nomials approximating uniformly @ @ on each compact set 


of C and Gin ~ ~ ce | Gin) |» therefore the result. An imme- 
diate application is the weak central limit at order n. 

Lemma 2.18: Z,=p sO Kei ) converges 
weakly (p— + œ ) for the functions of H„ (C). Thus, for 
all functions geH,,,(C), one has 


x 1 = 
EPZ) =z; f P(A )exp a 
+ go) 
Proof: 
Flexo 218 
(Y eC), Efexp(AZ,)}= II Bex» ) 
k=1 
Pp A” 
sil ett 
=exp(å"/n!). 


Then, by using the results of Lemma 2.17, one obtains the 
first identity. For the second identity, it is sufficient to 
develop the analytic expression of @ and to conclude by 
remarking that E{Z>7*} =6, of (ng)!/q!(m!)4, for all geN 
and r=0,...,.2—1. One remarks that E{o(Z, )} is inde- 
pendent of p, i.e., it is a stationary sequence if it verifies 
that U,=ZR_op {0} /[k(n!)"] is convergent (p 
++). 

If peH, (C), U, is convergent. If pH, (C), one can 
give an extension to this result by using a singular integral 
to give a formal sense to 


1 
x | PlA)exp 


Corollary 2.19: (Y @H,,(C)); one obtains 


(—iA)” 
ni 


Jaa 


E{9(Z,)} =E{P( Gin) } = EL Pol | Gin |)}, (29) 


with G(x) =(1/n)(2%2)9(x@,)). Noting that Z=w 
—lim Z, One obtains 


p> + o 

(Y AeC), Efexp(AZ)}=exp(A"/n!), (30a) 

E{|Z|?}=+ 0, (30b) 
Law 

Z, ~ Yo+Ginp (30c) 


where Y, and Gin are independent rv, Y, p being C valued 
and verifying 


w— lim Y,/p "72/0 


P>t+t o 


in the sense of tight convergence”! (G, and G, being two 
independent real standard Gaussian rv). 

Proof (see Ref. 13 for a complete proof): Here (29) 
results directly from Lemma 2.18. For (30), the idea of 
the proof consists in the study of the sequence of rv 
Z,/pr—7"/O") fie, = Re(Z,)/p'"—-/2" and 
Im(Z,)/p" 7? C®], and then, to “conchide via he clas 
sical central limit theorem. 

Corollary 2.20: More generally, let (X;) n+ be a se- 
quence of independent real positive rv with the same dis- 
tribution, verifying 


O<E{Xi} and ELX* <+ o. 


Then, if ( ef} ) iene denotes a sequence of independent 
Rademacher rv of order n and independent of the X, one 
has T p=Pp 7" KONA 161} X; ;), which converges weakly (p 
++), in the previous sense to T, with T verifying 


Law Vn je 
T ~E "Xi }Gin tY, 


where Gj, is a standard n-Gaussian rv independent of Y, 
and Y is a C-valued rv satisfying 


(Y keN*), E{Y*}=E{¥*}=0, 
E{|¥|7}=+0, 
Law 


Y~o,Y, for all g=0,.....—1. 


Proof: That is an extension of the classical central 
limit theorem. In the case of signed or complex measures, 
the result has already been proved by Refs. 9, 22, and 2. 
Some expansions are viewed in Ref. 13 and, in particular, 
the study of the Sring series of order n, i.e., of the 
form S,= 2h. 1pefa] , where (@k)ken+ iS a Sequence of 
l, (R T and Gi )kenx is a sequence of independent 
Rademacher rv of order n. 
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Now, one shows an interesting asymptotic behavior 
for the |G,;|’s law, which is given by the following. 
Lemma 2.21: (3 B„>0), such that 


(Y x>0), P{|G,,)| >x}<exp{—B,x"}; 


r—l1 


(n~!+n =1). 


Proof! Let Y,=n/[(n—1)]”"-). One has 


w 
|Gin| ~ Mi hnj (see Ref. 13), where (Hp ocjen—2i8 a 
family of independent real positive rv defined by 


(Y x>0), 


— u” 


n x, 
Phryx) = EDAG 1) /n) Í, wex ay 


n 


du. 


One verifies easily!> that all the moments of | Gin] and 
bln jare identical, and also the characteristics func- 
tions and the laws. 

Noting that C,j=n/(y/*?“"T((j+1)/n)), one de- 
duces that 


—y" 
du 


n 


n—2 
+o . 
P{|Gin | >x} = Il Cu | wexp 
j=0 xj 


n—2 +0 
<(m—1)! II ae 
j=0 xj 


n 


xexp du, 


n 


with M ixj=x, because the following functions F; are 
positive on R_: 


+00 ? 
F= f [Cn 1u"! — Cpp ]exp{—u"/y,}du. 
Indeed, 


+o , 
F(x) = f [Cnn 14"! —C, wJexp{—u”/y,}du 
0 


-Í [Cnn-14"7!— nj’ lexp{—u"/yYn}du. 
0 


From that, one deduces F;(x)>0 for O<x 
<lCnn-1/Cnj) OIP and = Fi(x)<O for x 
S (Canai Cr) TR one has, equally, F,(0) =O and 
F(x) 30 for x>1; thus Fj(x)>0 for all x>0. 


Finally, one gets 


n2 y =x} 
P{|G n—1)! — ex 
{| Gtny| >x}<( »| TH z= e( | 
j=0%} 


<exp =: 


With the help of the concavity of the Log function, one 
gets 


(V meN*) (V w,...,0,ER*), 


1 
Log{ (w+: +Wm)/m}>— {Log w+’: +Log Wm}. 


Hence one has P{|G,,)| >x}<exp{[—(n—1)/y,] 
xx”“—1)}, Detailed results are given in Ref. 18 on den- 
sities similar to p,(x). 

Remark 2.22: The studied sums of rv are not infi- 
nitely divisible in the Doob sense,’ because Gin is a RY". 
valued rv, whereas the sums of rv are C valued. Indeed, 
one has studied convolutions of measures defined on R1/”, 
which are constructed on (RY")?. Furthermore, all the 
given results are derived from the possible construction of 
the n—R_, -isotropic part ġo of an analytic function in C, 
¢, defined by 


1 n—l 
w=: 2 s(o)), (xeR). 


=0 


Formally, we have used results on the lacunary series of 
order n, i.e., SO (z) =Z oan" from the study of the 
entire series 


p 

S (0z) = 2 a,(@,z)*, for r=0,.„n—1. 
In Ref. 13, weak convergence is shown for S, on the 
C—C) continuous and bounded functions, by using an 
approximation on /,(R'”"). 


lil. THE STANDARD N-BROWNIAN MOTION BASED 
ON RY? 


One gives a definition of this process allowing his 
construction, following Levy’s ideas,*”> generalized at or- 
der n. 

Definition 3.1: Given a complete space of probability 
(0,F,P), where F is the ø algebra @ [B(R'/")]*+, with a 
filtration (F, t>0), one calls standard n-Brownian mo- 
tion based on Rye F,, adapted, all the stochastic processes 
noted down X/,;={(X,,)(¢);t>0):R,—+C} defined on 
[CU{ co }] F+, by the following assertions: 
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Xin (0)=0, as (32a) 
(V O<s<t<o), 
X (nyt) —X {ny (s) (32b) 


is a C-valued rv independent of F, (in fact, X,,) has inde- 
pendent increments); 


(Y O<s<t<+o), 


Xim(t)—X p(s) (t—s)!"Z, (320) 


Xin] has stationary and self-similar aos) ; 


(k) 


Z=w-— lim p [a] 


Pp +00 


) (see Sec. II); 


P 
al x G 
k=1 


w— lim [Xin (2) —X in (s))/(t~—s)""] 


tos 


exists as an n-Gaussian rv, (32d) 


[the differential increments of Xj, are R}/” valued and of 
order (dt)'/"]. 

Remark 3.2: The no-respect of the condition (32d) 
would imply a no-isonomy of all the processes, verifying 
the other assertions (see the structure of the o algebra). 
Furthermore, (32d) gives the fractal nature of this pro- 
cess that is totally different from the fractional Brownian 
motion of Lévy-Mandelbrot, © because his disjoint incre- 
ments are independent and his variance is infinite. Now 
one gives a construction of X/,) by a linear cut approxi- 
mation (time discretization). 


i Ml 


kT kT 
pA A A 
q 
kz 
X Gia] —I|>€], 
q 


P{|Z,(t) —Z,(s) | >)=P| F 
noting that k; =[qt/T] and k,=[gs/T]. Then, 


ky-1 


e 


r=k, 


l/n 
P{|Z,(t) —Z,(s)| >d=P| (7) 


e| 5, le) 
T\* /T 


7 


q 


kT 


q 


3135 


Lemma 3.3: Let (Z (t); te[0;T genx (0<T < +) 
be the sequence of linear cut processes defined from the 
uniform partition of [0;7] {kT/q}; k=0,...,q—1, and the 
family of independent standard n-Gaussian rv {Gin 
X (7/9) Frvqcto,11n@» by the following assertions: 


k-1 


—1/n 
alga) (3G) 
for all k=0,...,g—1. (33a) 
If te[kT/q;(kK+1)T/g[, then 
2,(1)=2,(—| poe IZ Ake T/A 
q (T/q) 
—Z,(kT/q)). (33b) 


Then, w—lim,.. + Z(t) exists in the sense of Lemma 
2.14, and has the same law as 7!/"Z [see (32c) for the 
definition of Z]. 

Proof: From his construction (Y O<s<t< +œ), w 
—lim,_, + wlZg(t) — Z,{s)] exists in the previous weak 
sense of Lemma 2.14 (see Ref. 13 for a detailed proof), 
and Z, has independent separated increments. Next, one 
will show that the rv Z,(t) (qEN*) are tight for a fixed 
te[0;7] (relative compacity in the sense of the Prohorov 
theorem?'). First, one needs to know the behavior of 
P{|Z,(t) —Z,(s) | >e}. 


Lemma 3.4: (Y gEN*) (WV €6>0), 


P{ Sup |Z,(t)—Z,(s) | >e}<exp{—B,67 /"~ Per’}. 


|t—s] <5 (34) 
Proof: 
T 1/n—1 ky kT T V/n-1 
(a) a a) 
qt kı qs ky\\ (T\~' 
T=) afa) -7 oim(Z)|(3) [> 


(E 
+) 


kz 


q 


qs 
T 


Ih 


Ce)l(F) 


Gin 
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Thus, 


P{|Z,(t)—Z,(s)| >€} 

kyl P 
Ht {lm (5) 
a (EG) Lom a)l] 
xa((F-«) (4) eala) 


with Bj+P,+2,l,a,=€(T/q)—¥™. 
By using the asymptotic estimation (31) on P{ | Gin] 
> x}, one gets 


>a, 


P{|Z,(t)—Z,(s)| >€} 


<! II exp(—B,07/"" n] 
jek, 
B B,(T/q) a 
G a 
B,(T/q) n/(n—1) 
xexp -B TT) $ 


That is a problem of variations with constraints, and one 
must minimize XZ_;z% under the constraint 22_a,z,=x 
(n’=n/(n—1)). For that, one notes that 


P P w 
Pee ae 


k=2 


the minima of F are reached for ðF/ðz;=0, i=2.,...,p; 
then the linear system in z,, 


p i 1/(n'—1) 
> QkZk+ (a (z) Jams, i= 1,2,....p, 
k=2 a; 

kÆl 


has a unique solution z;j=xa}"!/ (E-104), j=1,...p 
(the determinant of the system is different from zero, and 
one gets the coordinates z; of a certain gravity center): 
One deduces that 


E eid a a E ah 
È gxr — =»"( > at) 
k=1 (Bhat ay) k=1 


(this proof uses the same ideas as this one on Minkowski 
inequalities). Hence 


P{|Z,(t) —Z,(s) | > 


<exp{—B,(t—s)—/"—- Ye}, 


Corollary 3.5: One has that — lims_o+ 
P{Sup,,_5|<s|Z,(¢) —Z,(s)| >€}=0, uniformly with q 
for all €>0, i.e., the rv Z,(t) are tight and converge 
weakly in the sense of H,,(C), to an n-Brownian standard 
motion based on R'/”, Xin} 

Proof: From his construction, X|,) verifies (32a), be- 
cause Z,(0)=0, a.s for all gEN*. Here 


Z,(t)—Z,(s)= a k ae kı 
o-z] a) 
qs T 1/n—1 k 
+(F-4)(3) Sola) 
aNd kg! > 
+(3) [3 eela) 


This is a sum of independent standard n-Gaussian rv, 
and one gets 


t/T—k;)]” 
EZD- a || 


(T/q) 


daza 
(T/q) q) 

Since gt/T—k,<1, qs/T—k, <1, 
—k,)/(q/T) =t—s, one deduces that 


and lim (ki 


q= +o 


lim E{(Z,(t)—Z,(s))"}=t—s. 
q> +0 


Via the results obtained in Sec. II, one deduces that 


Law a 
Xin (t) —X nls) ~ (4-8) "Z. 


Since the family {Gj,4(7/q) }ocrcq With a fixed gq is a 
family of independent standard n-Gaussian rv, and from 
the construction of Z,, one deduces that for q great 
enough, Z,(t) —Z,(s) is independent of Z(t’) —Z,(s') 
for O<s’<t' <s<t<T; then, going to the limit 
(q> + œ), one deduces that X;,) has separated indepen- 
dent increments. For (32d), one has an inversion limit 
problem, and one must study the ratio {[Z,(¢) —Z,(s)]/ 
(t—s)'“"}; also, one comes down to the zidy of 
E{[Z 4) —Z,(s)]"/(t— s)}, because E!/"{(-)”} seems to 
be the adapted norm for the study of rv of type S, 
= aG 

Corollary 3. 6: Note that Xip (t) =w— lim, n4 o Z,(t) 
(in the sense defined previously); then w 
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— lim, + [EX (t) — Ximls)]/(t — s)'/"] is a standard 
n-Gaussian rv, and also, Rin valued [in the sense of the 
convergence on H,,,(C)], 

Proof: 


g [EO O _ 5) 
usA [n] q 


<E"” [Eae zN 


(t—s)” 


„ [[Z0 —Zq(s) k\]” 
aee | m E] 


Also, 


—X k,\ 1" 
pe fi, ) 


(t—s)'/" 
or, EV Ze} + EM" ae 
Va-1 : 
A ol- 
Then, finally, 
on feo 
(t—s)'”" m |7 


<E + KSL = \G 


(t—k,T/q) (=) Mam) r| 
Because of the construction of Xp} limg- + €; = 0 and 
for q large enough, one can choose ts, verifying t—s 
<T/q and t—s~T/q (q-+ œ). Therefore 


=1; 
q 


i (t—k,T/q) (T\'"-} 
asi (¢—s) ( 


then, for g=[T/(t—s)], one has 


Xim (t)—¥; (8) k 
l/n in] [n] E *2 
g | (sy om (7) 


n 


<A(t—s) 


[where A is a positive continuous function in 0, verifying 
A(0) =O]. 

Remark 3.7: To follow the Maruyama notation, !2++4 
one writes, in a symbolic way, X mD =S Gin (5) (ds)!” 
[where {Gin (s) }o<s<r) is a family of independent stan- 
dard n-Gaussian rv]. 

To obtain other results, one can consult Ref. 13. 
Next, one gives a construction of X;,) using the Ciesielski 
process,” and also, the Schauder functions. 


Lemma 3.8: Let {Gin (k271);g E N*;k e€ N*;k < 24 
be a countable family of independent standard n- 
Gaussian rv, (independent only for different values of 
k2~*). Then, one defines the following sequence of C- 
valued linear cut processes X = a(t):0<t<1) by 
(G) X(t) =tGiy (1); (ii) X; is linear in each interval 
[(kK—1)2~%k2~%] and continuous in ¢ for all we; 
(iii) X,41(2k2~ +) =X,(k2-%); and we see 


(iv) Xq41((2K—1)2- 4+) 
= A(2k—1)2~- 4+) 


$2-C+9 ME ((2K—1) 2-44 V9, 
We have the following results: 


P{ lim X(t) exists for all O<r<l, 
q> +% 


uniformly in f}=1. (35) 


Proof: 


P{ Max |.X,41(t)—X,(t) | >27 2} 
te[0;1] 


=P{ Sup | Gj_)((2k—1)2- 0+0) 2-C +9” 
1<¢k<2? 


>2-V(2n)} 


Thanks to the asymptotic estimation on P{[G,,)| 
>x} (31), one gets 


P{ Sup |Xq4.1(t) —X,(t) | >27} 
te[0;1] 
=P{ Sup |Gin((2k—1)27 +D) p20 + VAM) 
1<k<2? 


<2°P{ | Gin | >20 +420) 


<21 exp{— (n— 1 Jyp 2O HV ODJ- D), 
and a fortiori, (3 qoEN*), such that (Y q>qo), 


P{ Sup | X41 (4) —X (t) | >27 OC" }.< 272-4 = 2-4, 
te[0;1] 


comparing the growths at the infinity between the func- 
tions x and exp(—x*) (a>0; x++0). Since, 
Zgen2~?=2 and Eycn2~”") is convergent, also bounded 
by a certain gg € N*, one deduces, via the Borel—Cantelli 
lemma, that for all ¢>Sup(qo,9), 


P{ lim ( Sup |X;,(t)—X,(t) |) >2-” go} <A,274, 
k++ t€[0;1] 
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where A, is a positive real constant depending only on n. 
One remarks that the maxima of the processes X, are 
growing with q and are attained at the top. In particular, 


(Ve>0), lim P{for k>g and 
q> +o 


O<t<1;] X(t) —X (t) | >e} =; 


therefore the result and, ipso facto, one deduces the a.s 
continuity of the paths for the limit process ¥ in: 
Corollary 3.9: The X X jni ’s paths are a.s continuous. 
Lemma 3.10: Let nn 2—4(t) be the Schauder functions 
defined by 


Alt] (t)=0, for |t—k2~4|>2~%, (36a) 
AMR) (k279) =2~ Gt DM, (36b) 
heL, is linear in two adjacent intervals. (36c) 


Then, the following series is convergent to an n-Brownian 
motion Xi] in the sense defined previously, and one can 
write (V te[0;1]), 


ay AMZ] (t)G{q\(K279). 
(koda) 

(37) 
= lim X,(t) as. 


q> +o 


Proof: The independence of the separated increments 
is obtained by remarking that 


X,(t)—X,(s) 
=X, (1) —X,([2%]27%) 


+X{( [24s] +1)2-9—X,(s) 


[2%] 


+ S 


m= {29s} +2 


[X,(m2~?) —Xf(m— 1)279], 


for O<s<t<l. 


For the construction of X(t) —X,(s), one uses only 
distinct Gj,(k2~%) of these used for the construction of 
Xt’) —X,(s') (O<s’ <t'<s<t<1); also, via the inde- 
pendence of the G,,,(k2~%), one gets the independence 
for separated increments, and the result is valid, going to 
the limit (q— + œ). One has 


[2%] 


> Giny (7279) 
r[2%s] +2 


X(t) —X,(s)=2-/" | 
+(t— [2%2]2—9) G,.(( [2%] +1)2-9 


+(([2%s] +1)27—s)G,n) 
<(([2%s] +929]. 


Also, X,(t)—X,,(s) is of type 2,a,G{4) (a,>0). Via an 
extension of the isometry with the cone /,(R..), one can 
find the limit law of X,(t) —X,(s) (q> + œ) in the weak 
sense, previously defined in Sec. II. One has 


[2%] 
E{[X,(2) -x (91-24 ( 2 1) 
r= [2%] +2 


+ (t—[2%]27%) 


+(([2%s] + 1)27-4—s)n 


-t—s, (q>+o) 


and therefore 


Tar 1/n 
Xin (t) —X in (s) ~ 0—58) "Z 


(Z defined in Sec. II as the limit rv obtained in the central 
limit theorem at order n). But, one has the following 
“strange” properties. 


Lemma 3.11: 
lim {|X -X)| }=+% (O<s<t<+o), 
OES (38a) 
(WAeC), 


Eexp(aLX tn (4) —X nj (5) I} exp{a"(t—s)/nl}. 
(38b) 


E{|X in (t) —X tay (s) P}=+ 0, (38c) 


(V geH,(C)), 
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ELHAX ini (t) —X in (s))} 
to gink) 
k! 


(0)(t—s)* 


| du (38d) 


1 A 
aie f o(u)exp 


{where (x) =Sp exp(iux)d(x)dx]. 
Proof: 


290 -2/nT (2) 
E{|X¥;t)—X;(8) ~rr ES) >+, 


(q> +œ). 


The other results are directly validated through the 
previous constructions of X;„ and the properties of sums 
of n-Gaussian rv studied in Sec. II (i.e., the central limit 
theorem at order 7). 

Remark 3.12: Equation (38d) joins together X;„} and 


the Hochberg’s process X bal ? (that is defined only for 


p 
n=2p for convergence reasons), and that verifies the fol- 
lowing. 


For all functions ġy (R-R) in the Schwartz class, 


Eldar (t) xE (s))} 


1 A 
Sr i o(u)exp 


Lemma 3.13: Each of the following processes is a 
standard n-Brownian motion based on Ri’, 


2, 


uP 
Tp) (t—s)} du. (39) 


Xin (t)—Xia (s), 4820 (Additivity property), 
(40a) 


OX ini(t), t20, keN (n-R, isotropy), (40b) 


CX in (t/c"), t20, c>0 (n-Brownian scaling), 


(40c) 
PX in (1/t), t>0 (n inversion), (40d) 
orh "IX n (SEAL) —X ny (8) ], 
keN, h,t,s>0 (self-similarity). (40e) 


Proof: Equation (40a) is evident; Eq. (40b) results 
directly from the property established on Z in Sec. II. For 
Eqs. (40c) and (40d), an easy proof is given in Ref. 13; 
Eq. (40e) is the result of (40a), (40b), and (40c). 

Remark 3.14: In Ref. 13, one shows the strong Mar- 
kov property and the no-differentiability of the sample 
continuous paths of Xin» with the help of an evident gen- 


eralization of the case n=2.”> Next, one studies the vari- 
ations of order n of the a.s continuous sample paths of 
Xi, Which is the starting point for the definition of a 
stochastic integral and an It0-Taylor lemma of order n, 
using an extension of the isometry on the cone /,(R +) 
Lemma 3.15: Let T>0 and {t,,} be a growing par- 
tition of [0;7], more and more thin (i.e., O= go<*** < Egk 


<` <tyg=T; limy- + SuPock<ql é k+l T tak| = 0); 
then, one has 
q—! 

lim [Ag Xin] (41) 
q=>+o k=0 


exists and is R, valued [in the sense of quadratic mean 
and probability convergence; Ng kX in =X ny (gk 1) 
—X in (tg,x)]- 

Proof: Indeed, one has chosen X, [4] among the 
processes, verifying w—lim,_,,+ [[Xpqy(t) — Xims) (t 
— s)'/"Jis R'/” valued; then, if the limit of 3{Z)[A,, Xq]” 
exists, it is KR, valued, and to prove 
the quadratic mean convergence, it is sufficient to 
calculate E{((Z{Z[A,,-X,q]") — TY} in place of 
E{|(z fa0l Ag X wl”) —T| *}; therefore 


ol((Zeanr)-) 


k=0 


q-1 4 
=z|( p> (14g nl"—Aae!} | 


q—1 
= 2 ECLA gX tn ]”— Agt), 
because T=5{Z0 A git [Ag-Xty]"—Az,at is independent 


of [A, Xtal" — 44 for kÆj and Ex[Ag Xp)” — Agt} = 0; 
furthermore, 


E{( [Ag cX tn] }"—Agit) }= EL [Ag AX in) }*} 
—E{ [Axe] Y; 


then, one gets 


aC E eea) 


k=0 


g—1 
= È ATM g Xn PE Ag at?) 


(2n)! a 
act ( 3, (vt) 
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Qn) 4] 7C Sup [Acst]) +0 
<| can? 1] 7 Sup [Maat 
(q> +o). 


Corollary 3.16: (Y qEN*), 


A(x) 


k=0 
(2n)! S 
= son e [4,4] >0. (42) 


Proof? [(2n—1)(2n—2)-+:(n+1)V[(n—1) (n—2) 
*++2X1}]>1 for n>2; then the result is validated. 

Via the Markov—Tchebychev inequality, one deduces 
the convergence in probability for >> ia ae @ tal” 
(q+), and for a dyadic partition, one has the a.s 
convergence to T (in fact, that is the law of great num- 
bers). 

Corollary 3.17: More generally, for a partition 
{ta rto<k<q Of [s;¢] as the previous one defined in Lemma 
3.15, one has the following energy identity at order n, i.e., 


q—! 
EXMO Xm) È EAX nl") 
Š (43) 


Proof: The proof uses the independence of the sepa- 
rated Xis increments and the values of the different 
moments of A, ,X;, up to the order n, i.e., 


E{ [ Ag KX [n] 3 =6 jn [ Agt] , for j=1,2,..ħ. 
Notation 3.18: Formally, one notes that 


1 


J, 2Xn(¥=w= im E AX 44) 


g 
q>+o k=0 


(in the sense of the strictly definite positive functional 
E{C)}, with 7 = inf,ewe{np/j € N*}. 

One then deduces, with the help of the stochastic 
integral defined independently in Sec. IV, the following 
recurrence relations. 

Lemma 3.19: One has 


t i ; j=} t 3 
f [dX,\() nO- 1 CG f Xis) 
0 r=1 0 
X [dX in (9) 1’. (45) 
Proof: This identity is true, even going to the limit 


(g—> +œ) on the partition {t);}o<kc@ and it uses only 
the binomial identity. Furthermore, So[dX,,;(s)) defines 


a complex-valued stochastic process with differential in- 
crements as R?” valued (RY/" = {zeC;z=pa xs PER; 
O<k<n—1}). 

Corollary 3.20; Noting that po=rj/n and r 
= infpens{np/j € N*}, one gets the following. 


Ie pel, if [dX_)(s)¥=0 in the sense of E{(-)’}. 
(46) 


Proof: Indeed, [A Xin ~ [Ag at G 4 (t,,4) 
(q— + œ), and therefore the differential increments are 
RŽ” valued. Then, if po41, one has 


q—! r 
E ( > [Apain l’) | 
k=0 
q—-1 


my 2 Ma Xi Agr] 
beo K= 


Therefore 


q—!1 r g-! 
zl|( bs lapain) = $, E{[AgeXtny)"} 
k=0 k=0 


ye! 
ea ( > [Apat] >0, 


~ pol (nt)?o \ Ao 
and going to the limit (q—> + œ ), one has 


1 


E [Ag X in] r) | 


z 
( > 

k=0 
á (npo)! 


— tS A. ,t]20-!) 30 = f 
‘Do! (n!) SA qkt ~") (q> + 0) 


Corollary 3.21: In the case pp=1, i.e., n/jEN*, 


N(tq)—1 ; t ; 
wo tim È AX S OY (47) 
q>+œ k=0 0 


is a standard n/j-Brownian motion, [N (ż,q) is defined by 
tanap St < laNg 

Proof: One studies the sums of rv of type Z(t) 
= Sh AX inl’, where {t}, k}o<k <q is a growing par- 
tition of [0;7] more and more thin. 

Because of the independence of the X/,,’s increments, 
one gets easily 


(V O<s<t<T) (VY m=1,...,r—1) 


ELZ? (t) —Z (s)]"}=0, 


and by construction, Z‘? has independent increments, 
and therefore 
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ELIZ (t) Z? (s) 1} = taneg — taNas 72 —-S 


(g> +œ). 


Furthermore, 


Mop Ze? = [Ag X n ~ Aget] Ging) (9 +00). 


Using the previous results on the central limit theo- 
rem at order r given in Sec. II and the construction of the 
standard r-Brownian motion based on RY” , one deduces 
that Z!” converges weakly to a standard r-Brownian mo- 
tion. One can see Ref. 13 for the details of the proof and 
the characteristic properties of Xin} 

Remark 3.22: One remarks that the process, 


Z(t) =X (XG (t))+iXAX (H), (t>0) 


verifies 


(V O<s<t<+o) (V AeC), 
Efexp(A[Z(t) —Z(s)])} 
=Ey {Ey x, {exp ALX (Xo ()—X (X (s)) 


HiX (H) -XAA (SDH 


A? 
= Ey | exp (Fo -x80 


A? 
Xexp (-F 1x6 x5 001)] 


At 
zE] 
(48) 


[where Xo, X,, and X, are classical standard and inde- 
pendent real Brownian motions, X¥¢ (t) =Sup(Xo(2),0), 
and Xo (t)=—Inf(Xo(t),0), Ey being the average over 


A2 
exp (FPD —Xo(s) 1) 


=Ey, =exp 


But, the increment Z(t)—Z(s) is dependent of 
{Z(u); O<u<s}. To a multiplicative constant less, the 
increments Z(t)—Z(s) have the same extension of Fou- 
rier transform as the increments of YX, (t) —Xj4)(5), i.e., 
E{exp(A[Z(t) —Z(s)])} =exp{ (A*/8) (t—s)} and 
Efexp(A[Xjay(t) —X4(s)])} =exp{(A4/24) (t—s)}. 

The E X pX partial moments are as unbounded, con- 
sidering Xt (t) and Xo (t) as random times related, re- 
spectively, to X, and X,; Z is a fractional Brownian mo- 


tion of order 1/n,!? and one can say that it is constructed 
from the standard real Brownian motion by a stochastic 
iteration. 

Conclusion 3.23: Finally, one remarks that X;,) has 
been constructed from the definition of a discrete n mar- 
tingale,'? i.e., from the following process. 

Definition 3.24: Let (Q,F,P) be a probability space 
with a filtration (F;,) xen; then, the C-valued discrete sto- 
chastic process, M = (Mj) xen, Fz adapted, is a discrete 
n martingale, if and only if 


M(0)=0, a.s, (49a) 
(V keN), |E{Mi}| <+, (49b) 
(V keN) [Mk —M;]” is a real positive rv, 

(49c) 
(Y O<j<k) (Y p=i,...,.2—-1), 
E{ (M; —M;)/F;}=0. (49d) 


IV. A STOCHASTIC INTEGRAL AND AN 
ITO-TAYLOR LEMMA OF ORDER N 


Preliminary 4.1: One gives an elementary construc- 
tion of the stochastic integral of order n in relation to 
Xin) by using an extension of the isometry between the 
cone of rv (it is a discrete n martingale), of type Z,a Gly 
[where (@;) zens € l (R) and {Ge} rene is a family of 
independent standard n-Gaussian rv], and the cone 
1,(R,) via the functional E'”"{(-)"} definite strictly pos- 
itive on this type of sums of rv; then, one deduces an 
Ito-Taylor lemma of order n. 

Definition 4.2: Let Ej ((0;T]) (0< T< + œ), denote 
the set of step functions ([0;7]XQ—-R,), F, adapted 
and a.s bounded; then, one defines for all functions 
fF (u,w)eE;,([0;7]), the following stochastic integral: 


t 
f f* (u,@)dX [nj (u) 


p-1 
= 2 J> (ty) [X in (tee) —X [ny (te 1, (50) 


[where {t,} is a partition of [0;7], attached to ft by: 
(Vo ueltgte il, f* (4,0) =f* (tp@)]. 

Lemma 4.3: For all j=1,..,.2—1 (Y O<s<t<T), 
one has 


t j 
E ( J. f* (10) 4X mu) ) Je] =0, 
S 


(5la) 
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E ( f “P*(weo)dX m) Jr 


=E | J (f* (u,o))" du jeJ : (51b) 


Proof (Y j=1,...n—1), 


E ( fP tuwdm) Jr 


pi j 
=E ( È S ino LX lta Xin Ct] Jr 


p-l 
= È ELL ft (too) V/F} E {Xin (tes) 


—X in] (tg) V/F,}=0, 


aay for j=1,...n—1, E{[Xpy(te41) —Xiny(t)¥/ 
Fyd = 0, ft (tpo) is F, „ adapted and 


E {[X in (te, 41) —X pny te) [Xn Cee) 


—Xinj (ti) PEs} = 


for all r=2,..../; m+ +m, =j and kj <k,<°*: <k, 
When j=n, it remains only that 


p-} 


È, EUS oo ]"/F,} E UX (tea) 


—X in) (tk) "/F 3 


p-t 


= 2 ELU (tpo) F} (tki — tk) 


z z| if (f+ (uo) )" du/F, 


Lemma 4.4: Noting that Min (10; T]) the set of ran- 
dom functions f*+(t@) ([0;7T]XQ-—R,), verifying 
ELSTU * (u,@))" du} < +0, for all 
f*(t,@)eMjp([0;7]), one defines Sof* (u, o)dX y(u), 
te[0;7] [approximating f* by an a.s growing sequence of 
functions g, 7 €E ([0;7])], F, adapted in the sense of the 
following convergence: 


T 
lim £ f [f* (4) -gi (to)]” dt} =0 
0 


prota 


Proof: Using the arguments known and developed at 
order 2,*7° one deduces the existence of the stochastic 


integral of order n by the extension of the isometry on 
L,(dP xX dt), which is given by the strictly definite posi- 
tive functional E{(-) ”}. 

Lemma 4.5: One defines S§Xfy(s)dXjq(s) for all 
te[0;7] and keN by a time discretization. Then, one goes 
to the limit on the time step in the sense of the functional 
E{(-)"} (for the convergence). 

Proof: S5Xfj(s)dX{,(s) is also defined by the exten- 
sion of an isometry first, considering the sums of rv of 
type 212 4 ta (Lg Ag Xin], where {t ;} is a growing par- 
tition of [0;], more and more thin, and remarking that 


e|( 3 2 EE | 
-F E{X in (tg IE {Ag im ]"/F, J 
(nk)! ( % 
k(n)" al. pee lagal) 


(nk)! t 
-aar (J, ut du) (q> +o). 


Corollary 4.6: One has 


t n 
E |( f° Xidm) | 


Corollary 4.7: Noting that Sy 
= DIG Xba (ty pA, ;Xin], one obtains 


(nk) ik! 
(k+ED In)’ 2) 


(Y @geN), E{(S541—S,)"}>0, (53a) 
(VY veEN), lim E{(So4.»—Sq)"}=0; 

q> +œ 
Sup E{S3} <+ o, (53b) 
gEN 
S, is convergent in the sense of E{(-)”}. (53c) 


Proof: In order to construct {t41 j} from {tj}, we 
put on a real 441 Joti Verifying 


t 


Qjo tat ljo Stat Lig t 1 < lgti = 49+ htr 


Note that t =t jy 25 tq+1 + 4=%y,41- Then, one has 
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So 1—Sg= —X ny (t) [X pny (43) —¥ ny (+X fy 
X [Xin (t2) -Xin (41) ] 
+f) Xin (43) — Xin (0) ] 

= [Xin (ts) —X pny (t2) XE (2) XC). 


=a LLMM 


Therefore, by using the binomial formula written under 
the form 


a‘ —b*=(b+a—b)*—b* 
k 

= > Œb" (ab) (a,beC), 
r=] 


one obtains 


EX (S341 -S I= ELX ny (t3) —X pny (to) "XP (22) —X ny (4) 13 


k n 
=z |( 2 CX Cti) [Ags 1X tay] ‘ [Ape Xin] 


k 
-( Dep cee rt DJELA gyin a) E{[Ag41,j41X[n]"}>0. 


Piep =l 


Then, using the Xis and AXj„?}s moment values, one 
obtains 


E{(S941—Sq)"} 


k 


= [Ag+ 1j+14] x pA Cp CR 
Cie 


[n(k—m)]}! 


XR mn mi 


(nm)! 
OPLA gr il” ) 


Finally, one carries out step by step, adding an interme- 
diary point at each iteration, and one gets 


v k 
EUS, -50"}=E|( pA ( 2 CIX ia Catwjo) 
X [Ags wX in] r) 


n 
O } | 


To develop this expression, one remarks that one ob- 
tains products depending of the parameter R<n of the 
following form: 


>—=—.-/$}])-—s XS 


i (> 


Z+°'+ZR=N 
1<w, <t < WRV 


Zi 
x [Ags wpigk [ny] ‘| [Agt wait Xin] 


C X (tat wrip) 


In such a product denoted by Pp, [A, q+Wp dot 1X [a]? is 
independent of the rest of Pp, for 1<zp<n; then PR=0 if 
Zpz=n, because 
E{ l Agr wp Xin] t= Sigal Agwat] 
for k=1,...,n. 


Therefore E{(S,,,—S ,)"} has the following simplified 
expression: 


E{(S94..—Sq)"} 
= » al( 3 2 CX (gt wip) [Ags wX [ny] r) 


x [Ag+ wit 1% tn] "| 


v 


= 2 [Apt mit itl 


w= 
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k 
i [n(k—m)]!(am)! 
x aaj ts Ce ‘Ce (k— m)!m!(n!) 
EEE 


x Bal Moemafl” | >0, 


for all g,veN*; furthermore, one gets 


E{ (S +o — Sa) "} 
<K( n( > [Avsnafl?] 


<K'(T)( Sup [A,#])-0 (g++). 
O<j<q 
Also, one has an order on the S, in the sense of the 
functional E{(-)"}, and since S, is bounded in the sense 
of E{(-)"}, one has a kind of dominated convergence 
theorem and S, converges in the sense of E{(-)"} [S, isa 
Cauchy sequence in the sense of E{(-)”}]. 

Corollary 4.7: From the existence of stochastic inte- 
grals of type Si ox ta(s)aX, tal(S), one deduces the existence 
of stochastic integrals of the following type: 
SAX in (S))dX in (s), where f verifies fand f” are analytic 
in C 


+0 
fz)= È az! (apze€), (54a) 
kt-2 
lim K o Onki) (54b) 
k++ nk 


|| f AK OX (| 


+? (nk)lan, f t\* 

È testi (a) | 

and the condition given in the corollary is sufficient to 
have the convergence of this series. The details of the 
calculus are given in Ref. 13. 

Remark 4.8: One can define a “local” stochastic in- 
tegral, i.e., $ a Tof(X tnj(s))dX jn (5) with the help of the 
stopping time T,=inf{T>0;|Xj4(7)|>a} (a>0), for 
an analytic function f in B(0,a)={zeC; |z| <a} (where 
f(z) =} Raz" is verified: the entire series S(z) 
=SEB[(nk)!a „a/(k+1)!](z/n!)* has a convergence ra- 
dius R>T (T< +œ)). 

Remark 4.9: Other properties of the stochastic inte- 
gral of order n are given in Ref. 13. Next, one examines 


the construction of the It0-Taylor lemma of order n, 
which will be written under a symbolic form (semigroup 
and random drift) with the following formula: 


AL FAX ta (0) 


d. 
exp (dX m OF) -t AX (OD 


n 


1 
= È Gh? EmO m] (55) 


(where fis in the previous class of analytic functions in C 
in a first approach). 

Lemma 4. 10: w—lim, .+ EIA, X inl is denoted 
by SoldX, in) (s)]* (with the previous conventions and no- 
tations, the convergence being defined in the sense of the 
functional E{(-)"?0”*}, po = infpen*{7p/k e N*}). 

If po=1, i.e., n/kEN*, then 


t 
i [dX in (s) ] =X in (t) (56a) 


is an n/k-Brownian motion. 
If po> 1, then 


f [dX in (s)]*=0. (56b) 


Proof: ït is sufficient to show the result just in the case 
k>n. One also has 


q—1 


qr k <t( Sup |AgXtnj|""") +0 
O<r <q 


(q> +o), 


since the sample paths of X;,) are a.s continuous and {f, ,} 
is a growing partition of [0;7] more and more thin. 
Lemma 4.11: For all meN*, one has 


—k 


1 
Xin (= 2 H 


(57) 


Proof: Using the binomial formula and the discreti- 
zation in time of the stochastic integral, one gets 


Xin (t r+) = [Xin (2, + Ag Xin)” 


=È CX (tar) Bg X tml 


therefore 
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Xin Cars.) -Xn lard = È CX te Ca Me Xml 


and summing on r, it comes to 


m 


Xin (D) = 2 al 3 Xia “(inl Aget ini)” > 


then going to the limit (¢— + œ), one obtains with the 
previous notations 


mide È ct S (s) [dX in (8)]* 


De 
Tik 


df k k 
Slx (891| dX im (0016 


Lemma 4.12: For all functions f analytic in C verify- 
ing the previous conditions, one has 


n 1 t k 
AX in (0-0) = 2 a J, fF (Kin (5) 


X [dX in (s) 1% (58a) 


i.e., with the differential notation 


d. 
dl FX ON = [exp (Xn OF) EmO 


n 


1 
È Bl Am OANA ml) 


(58b) 


Proof: The Itô-Taylor lemma has been proved for the 
polynomial functions of type P(x)=x™” (meN) in 
Lemma 4.11; therefore it is proved for the polynomial 
functions with complex coefficients. To deduce this 
lemma on some analytic functions, one uses the Stone- 
Weierstrass theorem of approximation of an holomorphic 
function in a compact set of C by a sequence of polyno- 
mial functions in the sense of the uniform convergence. 
First, let f be analytic on C and {Pyatmen (a>0), a 
sequence of polynomial functions approximating f on the 
compact set B(0,a): One has (Y zeB(0,a)), 


-1 


z” 
Pmal2)=f(0) +z O) rT O 


2 (z—2z')"- 1 cn 
+f “Gat fm Pmnalz az’. 


Via the Ite-Taylor lemma of order n on Pn,@ one obtains 


PrnAX in] (t))—Pra(0) 


Pm (s))dX inz (8)- 


Then, going to the limit (m— + œ), even if one might 
stop Xin] on B(0,a) with the help of the stopping time T, 
one gets 


SX in (tA Ta))— f (0) 


n 


1 taT, 
-2m lim f * POA MAX 


t tm +o 


Since 


lim E| 
m— + œ 


PRE D-OE O) Me as 


0, 


one deduces the following result on the stochastic inte- 
grals: 


AT, 
tim [° POX pDA pay (8) * 


m— +o 0 
IATa o 
= f fF} (Xin (AX inz (8), 


in the weak sense given previously, defining the stochastic 
integral for all a< +æ. Since lim, Ta™= + œ, one 
then goes to the limit (a—> + œ), on each side of the 
identity. Also, one obtains the Itô-Taylor lemma of order 
n for the analytic functions in C, verifying 


lz | ig LAPE)" | iver 


which gives the same type of existence condition as that 
obtained for the construction of the stochastic integral of 
order n of f™. 

Remark 4.13: In Ref. 13, one goes further, studying 
the diffusions of order n of type 


n 


dé(t)= È a,(t,a)dX\(t), 


i=1 


where {Xfi ti }i<icn is a family of independent i-Brownian 
motions. One therefore obtains an It6-Taylor lemma on 
these diffusions. 
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V. RESOLUTION OF HEAT PROBLEMS OF ORDER 
N; MEASURES CONSTRUCTED FROM R” 


Definition 5.1: The heat polynomials of order 
n,>627,28 W,,p(t,x) are defined by the generating function 
exp(xz+ tz"), i.e., 


(V x,z,teC), 


+o 2? 
exp(xz+iz")= > Wrpltsx)— , (59) 
p=0 P: 


and by an identification on the formal series, one obtains 


k 
xr 
W,(ox)=pli} > =A 
nol ISP s Ar 


= E{(x+ (ntt) "Gin YP} 
=E{(x+ (n) "X in (OP } (60) 


(where G,,) is a standard m-Gaussian rv and Xj,) is a 
standard n-Brownian motion). 

Corollary 5.2: W,,,(t,x) satisfies the following heat 
problem of order n: 


7t. 


ðu "u OTIxe 
Bp bb) = gatx), (t,x)e[0;7] XC, 


u(0,x)=x? (xeR). (61) 
Proof: 
rex) =p( wey Z5) 


with m=r+1 and q=k—n. Thus the result is proved. 
Corollary 5.3: More generally, let 


n k—1 
(—1) 

> k Akxy, 

k=1 


g,(A) =exp (A,x eC) 


and {H,(x),....X,) }yeN denote the sequence of polynomi- 
als with n variables defined by 


+ 00 


galà) = 2 Hí Xis-Xn)AP, 
P= 


i.e., 
1 (dg, 
H,=— (a)| 
P p (e ieo 
Then (Y peN), 


H, ( f dka. j CAEN f [AX tn (8) 1") 


t ti byl 
= fatale if AX yy (ta) °° f dX in (tp) 
(62) 


(in the sense of the stochastic integral of order n iterated 
at order p). 
Proof: It is sufficient to verify that 


n 


(=1)4-'! t 
W(t) =exp » ep f. Aa (914 


is the unique solution of the stochastic differential equa- 
tion of order n: 


dM, (t) =M, (t)dX in (8), 
M,(0)=1. 


Thus, writing the present formal series, one obtains 
the following recurrence: 


Z,(t) = f Zp—1(5)4X in] (s), 


where 


Z,(1) =H; ICTO f [4X 1m (8)]") 


For the details and an extension of this result, one can see 
Ref. 13. 

Remark 5.4: Hochberg? has already given the begin- 
nings of this result, and via the Itô-Taylor lemma, the 
H,(x},....%,) are explicitly defined. 

One can equally notice that Jumarie” already linked 
g,(A) to the following evolution PDE of order n: 


Ou x ak ou 4 c 
g= i xk box) (t,x) €[0;7] x 


k=1 (63) 


u(O,x)=f(x) (xeC) 


J. Math. Phys., Vol. 33, No. 9, September 1992 


Phillipe Sainty : Brownian motion of order n 3147 


(where f is a given function with “good” analytic prop- 
erties). 

He uses an entropy notion. To solve this problem, 
one constructs stochastic processes called semi- 
martingales of order n” in place of X [np defined by 


n 


aX(t)= È a AX {A} (0), 


with ¥(0)=0 a.s, a,eC and (XA Jicken a family of in- 
dependent k-Brownian motions. 

Lemma 5.5: Let f be a function verifying the condi- 
tions of the Itô-Taylor lemma of order n (Lemma 4.12). 
Then, one gets the following Dynkin’s formula: of order 
n, ie (Y a,EC), (V AteR*), 


Ef f(x +0,X im (t+ At) )}— ELS (SHAX int 


1 t+At 


=- apE{ f(x+-0,X in (8))}ds. (64) 
t 


n! 


Proof: The proof given in Ref. 13 uses the Itô-Taylor 
lemma of order n (Lemma 4.12), i.e., 


Sx pX in (t+ At))—f(e+anX pn (2) 


tA 


ATN! 
=È ñj, © af (x + ară n (SAX n (8) 1%. 


Then, going to the average on each side of the identity, it 
comes out that 


EX f(x +-0,X in (t+ At))—S+0,X in (t))} 


] t+Ar 
Thn! f 


since all the other averages are null for all xeC, t>0, 
AteR*. 
Lemma 5.6: The heat equation of order n, 


an EL f(x +aX[n}(s))}ds, 


nou 
etx) = Tl Ox saltx), (tx)e[0;T] XC 
(65) 


u(0,x)=f(x) (xeC;a,eC*) 


(where fis a (R-R) function, continuous in L,(R,dx) 
with a holomorphic extension to C defined by f(z) 
=1/27Spf(A)exp( —iAz)dA, with 


Fa) fi exp(iAx)f(x)dx (AeR)), 


has a solution given by 


u(t,x) =E{ flxt+a,X in] (t))} 


t 
—idx-+(—iayA)" =} dA 


1 A 
=z f f(A)exp 
(66) 


Proof: Going to the limit (At+0) in Lemma 5.5, one 
deduces immediately the result. 

Remark 5.7: There is a relation between the class of 
admissible functions, a, n. For n even, one recognizes 
the Hochberg’s result,” and one can choose a, 
=exp{im(1/n—})}, i.e., až= —(—i)" to obtain the big- 
gest class of admissible functions. If one chooses a,, 
=exp{(ir/2)(1/n—1)}, ie. ati"), always for n 
even, it is sufficient to have f(x)eZ,(R,dx); also one can 
take fin L OL, by using the Plancherel identity. For n 
odd and a,=1, one can take equally fEL NL. One can 
equally remark that the formula (66) gives the solution 
with the help of a pseudodifferential operator (Parseval 
identity, I Ref. 13). 

Extension 5.8: Using the same method as in Lemma 
5.5 with a stopping time, one obtains a Dynkin formula of 
order n with a stopping time and in Lemma 5.6 one can 
extend the result to holomorphic functions fin B, (a>0) 
defined by B,={zeC; | Im z| <a}. By a more general way, 
one approaches fon a compact set of R by a sequence of 
polynomials {P,},cn in the sense of the uniform conver- 
gence, and the formal solution is given by 


u(t,x)= lim E{P(x+a,X [nj (t))}- (67) 


q= +o 


This solution is admissible for n even, CeO 1/n 
—})} and feL,MNL, and continuous in R. 


Appendix 5.9: In Ref. 13, one studies more generally 
the PDE of order n via the notion of stochastic differen- 
tial equations of order n, in the scalar and finite- 
dimensional case, which will be the subject of a next pa- 
per. Moreover, one will develop the notion of discrete and 
continuous Martingales of order n. Next, one gives some 
informations on basic measures implicitly used in this 
paper to define the standard n-Brownian motion. We will 
content ourselves to describe here the case n=3. 

Definition 5.10: Let v denote a probability measure on 
(R,B(R)), with a density v'(x) in relation with the Le- 
besgue measure defined on (R,B(R)). From v, one defines 
the measure pz on (C,B(R'?)) by the following assertions: 


BRI? ) is the o algebra generated by sets of the following 
type: 


i13(x,h) ={zeC;(z—x)*e[0;4°[}, where (x,h)ERXR*. 
(68a) 
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Let 1;={i1/3(x,h);(x,h) ER X R¥ }, 
and denote by [x;x+A?[ =(i1/3(x,4))°. (68b) 
Then, u verifies (Y A,Bel;), 

u(A)=v(4?), (69a) 
(ANB) =v( ANB). (69b) 
Lemma 5.11: Immediately, one deduces that 
(V A,Bel;), p(AUB)=v(42UB%). (70) 
Proof: 
(AUB) =p (A) +4 (B) —p(ANB) 

=v(A3) +¥(B?) —v(43 NB) 

=v(AUB). 
Remark 5.12: First, 


C= U U inlxan!®), 
neN* x€Q 


and then, by using (70), one has 


w(C)=v( U U [xgx,4+2[)=v(R)=1. (71) 
neN* xEQ 


Equally, one gets the following. 


Lemma 5.13: 
u(R)=1 and p(C/R)=0. (72) 
Proof: 


R=N U in (xpg ^) 
gen* X€Q 


because lim g~ 4 
qa +00 


+o "I [(k=1,2) have a length tending to 0 (q 
++), while R C Uxegit3(xsq7') for all geN*]. 
Then 


w(R)=p(N U iys(apq7")) 
gen* XER 


= lim u(U i13(x5q7'”)) 


q>+o xQ 


= lim v (U [xgx;+q7']) 


q> +o x 


= lim v(R)=1. 


q= +œ 


= 0, and then the branches [x;x; 


Then one obtains u (C/R) =0. 
Lemma: 5.14: More generally, let a,b€R (a <b); then 


H([a;6[) =v([a;6[). (73) 
Proof: One uses the previous method with 


[a;b[= N U 1/3), 


gen* xEQN [abl 


i173(X07 


Lemm a 5.15: (N FeB(R'/3)), 
BF) =0>v(FNR) =0. (74) 


Then, the parts of B(R'/*) that do not contain a no 
empty interval of R are negligible in the Lebesgue sense. 
Lemma 5.16: Let 


q-1 
J= U iysla+(k/q\(b—a), [(b—a)/g]'®), 
k=0 


with [(b—a)/q] <1 and qeN*. Then 


u(Ja)=v([a;b[), (75) 
and going to the limit (¢— + œ ), one will write formally 


du 
qay”! (76) 


[in the sense of the Radon—Nikodym measure derivative 
on B(R'/) for u and on B(R) for v]. 

Remark: 5.17: Via nonstandard analysis, all the re- 
sults developed here can be rewritten following the ideas 
of Keisler!! (he has constructed a nonstandard Itô cal- 
culus of order 2 with the help of an infinitesimal step time 
discretization). 


VI. CONCLUSION 


We have given some indications for the development 
of a stochastic analysis of order n (n>2) and their ap- 
plications to the study of evolution PDE of order n. The 
singular properties of the sets 4Æ(0,h) generating 
BCR”) conduct to a “strange” behavior, in comparison 
with the standard real Brownian motion of Levy. In a 
future paper, we will develop the notions taken up here in 
the more general context of martingales and diffusions of 
order n. Applications will be given to solve the master 
equation and others.?*?%! Finally, one remarks that Za- 
chary”? has implicitly used measures based on Ri" to 
solve the Gelfand—Levitan problem of order n. 


'Block, Ark. Math., Astron. Fys. 7, 13 (1912); 2, 7 (1912); 3, 1 
(1912); 4, 8 (1912). 

? Hochberg, Ann. Probab. 6, 433 (1978); Proc. Am. Math. Soc. 79, 298 
(1980). 

3Doob, Stochastic Processes (Wiley, New York, 1953). 


J. Math. Phys., Voi. 33, No. 9, September 1992 


Phillipe Sainty : Brownian motion of order n 3149 


*Lévy, Processus Stochastiques et mouvement Brownien (Gauthier- 
Villars, Paris, 1965). 

5Kemnitz, SIAM J. Math. Anal. 13, 640 (1982); 16, 379 (1985). 

6 Watzlawek, Appl. Anal. 18, 283 (1984); M.h. Math. 100, 67 (1985). 

Motoo, Stochastic Analysis and Applications, Advances in Probability, 
edited by Pinsky (Dekker, New York, 1984), Vol. 7. 

® Nishioka, Jpn. J. Math. 11, 59 (1985). 

*Gaveau and Sainty, Lett. Math. Physics 15, 345 (1988). 

10 Mandelbrot and Ness, SIAM Rev. 10, 422 (1968). 

l! Keisler, Memoirs of the AMS, No. 297, 1984, Vol. 48. 

"Berger and Sloan, LNM 932 (Springer-Verlag, 1982), Vol. 1; Mem- 
oirs of the AMS, No. 266, 1982; J. Funct. Anal. (1983); J. Diff. 
Equations 57, 15 (1985); 62, 252 (1986). 

13 Sainty, Ph.D. thesis, University of Paris, 1987. 

“Ungar, Am. Math. Month. 89, 688 (1982). 

SVan Hamme, Am. Math. Month. 1, 70 (1984). 

'SMittag-Leffler, Acta. Math. 29, 101 (1904). 

Montel, Leçons sur les Fonctions Univalentes (Gauthier-Villars, 
Paris, 1933). 

18 Zolotarev, Transactions of Mathematical Monographs (American 
Mathematical Society, Providence, 1986), Vol. 65. 

19 Pollard, Bulletin of the American Math. Society 1948, p. 1115. 


Feller, An Introduction to Probability Theory and its Applications 
(Wiley, New York, 1966), 2nd ed., Vols. 1 and 2. 

?! Billingsley, Probability and Measure (Wiley, New York 1979). 

22K rylov, Sov. Math. Dokl. 1, 760 (1960). 

3K night, Mathematical Surveys of the American Mechanical Society, 
No. 18, 1981. 

2% Chaturvedi and Gardiner, J. Stat. Phys. 17 429 (1977); 18, 501 
(1978). 

BI and McKean, Diffusion Processes and their Sample Paths 
(Springer-Verlag, New York, 1965), Band 125. 

?6Gihman and Skorohod, Stochastic Differential Equations (Springer- 
Verlag, New York 1972), Band 72. 

27Feinsilver, LNM 696 (Springer-Verlag, New York, 1978). 

28 Nevai, 2nd Conference on Approximation Theory (American Mathe- 
matical Society, Providence, 1980), Vol. 3. 

2 Jumarie, J. Math. Phys. 26, 1173 (1985). 

Van Kampen, Stochastic Processes in Physics and Chemistry (North- 
Holland, Amsterdam, 1981). 

31 Dubinskii, Russ. Math. Surveys 37, 5 (1982). 

32 Zachary, J. Math. Anal. Appl. 117, 449 (1986). 

3 Diaz and Means, Ann. Mat. 97, 113 (1973). 

Roetman, J. Diff. Equations 9, 335 (1971). 


J. Math. Phys., Vol. 33, No. 9, September 1992 


